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In this paper, the dynamical attractor and heterochnic orbit have been employed to make the 
late-time behaviors of the model insensitive to the initial condition and thus alleviate the fine-tuning 
problem in the torsion cosmology. The late-time de Sitter attractor indicates that torsion cosmology 
is an elegant scheme and the scalar torsion mode is an interesting geometric quantity for physics. 
The numerical solutions obtained by Nester et al. are not periodic solutions, but are quasi-periodic 
solutions near the focus for the coupled nonlinear equations. 
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I. INTRODUCTION 



The current observations, such as SNela (Supernovae 
type la), CMB (Cosmic Microwave Background) and 
large scale structure, converge on the fact that a spatially 
homogeneous and gravitationally repulsive energy com- 
ponent, referred as dark energy, accounts for about 70 % 
of the energy density of universe. Some heuristic mod- 
els that roughly describe the observable consequences of 
dark energy were proposed in recent years, a number of 
them stemming from a certain physics [l| and the others 
being purely phenomenological '2'|. About thirty years 
ago, the bouncing cosmological model with torsion was 
suggested in Ref. [![, but the torsion was imagined as 
playing role only at high densities in the early universe. 
Goenner et al. made a general survey of the torsion 
cosmology Q, in which the equations for all the PGT 
(Poincare Gauge Theory of gravity) cases were discussed 
although they only solved in detail a few particular cases. 
Recently some authors have begun to study torsion as a 
possible reason of the accelerating universe 0]. Nester 
and collaborators [1] consider an accounting for the ac- 
celerated universe in term of PGT: dynamic scalar tor- 
sion. With the usual assumptions of homogeneity and 
isotropy in cosmology, they find that torsion field could 
play a role of dark energy. This elegant model has only 
a few adjustable parameters, so scalar torsion may be 
easily falsified as " dark energy" . 

The fine-tuning problem should be one of the most im- 
portant issues for the cosmological models, and a good 
model should limit the fine-tuning as much as possible. 
The dynamical attractor of the cosmological system has 
been employed to make the late-time behaviors of the 
model insensitive to the initial condition of the field and 
thus alleviates the fine-tuning problem Q. In this pa- 
per, we study attractor and heteroclinic orbit in the tor- 
sion cosmology. We show that the late-time de Sitter 
behaviors cover a wide range of the parameters. This 
attractor indicates that torsion cosmology is an elegant 



scheme and the scalar torsion mode is an interesting geo- 
metric quantity for physics. Furthermore, there are only 
exact periodic solutions for the linearized system, which 
just correspond to the critical line (line of centers). The 
numerical solutions in Ref.Q are not periodic, but are 
quasi-periodic solutions near the focus for the coupled 
nonlinear equations. 



II. AUTONOMOUS EQUATIONS 

PGT P based on a Riemann-Cartan geometry, allows 
for dynamic torsion in addition to curvature. The affine 
connection of the Riemann-Cartan geometry is 

r^." = r^." + \{V + + T\^) , (1) 

where P^^" is the Levi-Civita connection and T^^ is the 
torsion tensor. Meanwhile, the Ricci curvature and scalar 
curvature can be written as 



R = R + 2V^T'^ + -(T^.«r^'''= + 2r^,,T«^'' - 4T^T'') 



(2) 
(3) 



where .R^j^ and R are the Riemannian Ricci curvature 
and scalar curvature, respectively, and V is the covariant 
derivative with the Levi-Civita connection and T,, = TY,,,. 
According as Ref. [0| we take the restricted form of torsion 
in this paper 



(4) 



therefore, the gravitational Lagrangian density for the 
scalar mode is (For a detailed discussion see Ref.@) 
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{T^a^T'^'"' + 2T,,^T^"'' - 4T^T^) , (5) 
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Since current observations favor a flat universe, we will 
work in the spatially flat Robertson- Walker metric. Ac- 
cording to the homogeneity and isotropy, the torsion 
should be only time dependent, so one can let Tt{t) = 
and the spatial parts vanish since we have taken the 
restricted form ((4|) of torsion. For the general form, the 
torsion tensor have two independent components 
From the field equations one can finally give the neces- 
sary equations for the matter-dominated era to integrate 
(For a detailed discussion see Ref.jl]) 
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flQ and the energy density of matter 
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One can scale the variables and the parameters as 



t - 



H -> 



IpHo^, R IpHgR, 



(10) 



where Hq is the present value of Hubble parameter and 
Ip = vSvrG is the Planck length. Under the transform 
pO)l . Eqs. ([6])- ([8]) remain unchanged. After transform, 
new variables t, i?, <i> and R, and new parameters Cq, 
ai, and b are all dimensionless. Furthermore, the New- 
tonian limit requires ao — —1. Obviously, Eqs. ©-(IHl) 
is an autonomous system, so we can use the qualitative 
method of ordinary differential equations. It is worth 
noting that in the analysis of critical points, Copeland et 
al. [lO] introduced the elegant compact variables which 
are defined from the Friedmann equation constraint, but 
in our case, the Friedmann equation can not be written 
as the ordinary form, so the compact variables are not 
convenient here. Therefore, we will analyze the system 
of Eqs.®-® using the variables iJ, <f> and R under the 
transform (jlOp. 



III. LATE TIME DE SITTER ATTRACTOR 

In the case of scalar torsion mode, the effective energy- 
momentum tensor can be represented as 



-i^,H^+-{R- 



^)(3i/-$)2-Ai^^(ll) 
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-MR-R)-Tt'], 



(12) 



and the off-diagonal terms vanish. The effective energy 
density 



Peff ^ p + Pt = P + Ttt, 



(13) 



which is deduced from general relativity, p^f f = pt = T^ 
is an effective pressure, and the effective equation of state 



P + Ttt 



(14) 



which is induced by the dynamic torsion. 

According to equations ©-([S]), we can obtain the crit- 
ical points and study the stability of these points. There 
are five critical points {He, $c, Rc) of the system as fol- 
lows 

(i) (0,0,0) 

(ii) ((5(i+^_^)yB+C,|^/B+C,^2(i+iii)) 

(iii) (_(S™_A)^/B+C7,-i^/:B+^7,-«™) 

(iv) ((Hi+^+A)v/B^,|yB^,-Hi+^) 

(v) (15) 



i?2 _ 3a,H^ . (9) where A 



_ ^a2(l + ai)3(l + 9ai) 



8ai(l + ai) 

3^af (l + ai)3(l + 9ai) 



B ^_(l±iiill5±£fii) and C = 



. Consider that the parameter b is as- 
sociated with the quadratic scalar curvature term R^ , so 
that b should be positive . Evidently, the critical points 
{Hc,^c, Rc) are not real values except (0, 0, 0) in the 
cases of parameters 6 > and oi > or — 1 < ai < —1/9. 

If we consider the linearized equations, then Eqs. 
dH]) are reduced to 

H^J^R, $ = R = -^<^>. (16) 

6ai 2ai b 

In the case ai > 0, there is only a critical point (0, 0, 0) 
for the nonlinear system. Eqs. (jl6p has an exact periodic 
solution 



-ff = ai?osinwt-|- -$ocoswt + _ffo - -*o, (17) 



P, 



$ = /3 ^ Ro sin cut + ^0 cos ujt, 
R = Rq cos u!t — /3$o sinujt, 



(18) 
(19) 



72ai ' V 

H{0), $0 = *(0) and Rq = R{0) are initial values. Ob- 
viously, (H, 0, 0) is a critical line for the linearized sys- 
tem. However, there is only a critical point (0, 0, 0) for 
the nonlinear system which is an asymptotically stable 
focus. In other words, there is no periodic solution for 
the nonlinear system since the corresponding eigenvalue 

is (0, ^ft*' \/ ^*)' Figdl we plot the orbits near 
the point (0,0,0) for the nonlinear systems. 

Next, the parameters are restricted within b > and 
— 1/9 < fli < or oi < —1 in this paper. To study 
the stability of the critical points {He ^c,Rc), we write 
the variables near the critical points in the form H = 
Hc + U, ^ = ^c + V, and R^ Rc + X with U, V, X the 
perturbations of the variables near the critical points. 
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FIG. 1: The phase diagrams of (//,$, J?) for nonhnear system 
©-(HI) with ai > 0. We take ai = 0.9,6 = 8 and the initial 
values (0.8, 0.5, 1.1). (0, 0, 0) is an asymptotically stable focus 
point. 

Substituting the expression into the system of equations 
©-(HI), one can obtain the corresponding eigenvalues of 
critical points (i)-(v) 

(ii) (-VB+C -3( '"+''i' -a) v'B+C-(i±|^ -3^)^/51^) 

(iii) ( v/B+C,3 ( Hi+iLLl ys+C, ( i±§^ -3 A) ^B+c) 

(iv) (-x/B^ -3( '"+''i' +a) ys^-(i±|ii+3A)\/B^) 

(v) (v'B3C,3(ai+iLa+A)^/B^.(i+|;ii+3A)yI^^ (20) 

The properties of the critical points are shown in tables 
U and ini We find that critical point (ii) is a late time de 
Sitter attractor in the case of —1/9 < ai < 0. 

IV. NUMERICAL ANALYSIS 

In previous sections, we have studied the phase space 
of a torsion cosmology. The de Sitter attractor indicates 
that torsion cosmology is an elegant scheme and 

the scalar torsion mode is an interesting geometric quan- 
tity for physics. In this section, we study their dynamical 
evolution numerically. The crossing of the w = —1 bar- 

TABLE I: The physical properties of critical points for ai < 
-1 



critical points 


property 




stability 


(i) 


focus 


±oo 


stable 


(ii) 


saddle 


-1 


unstable 


(iii) 


saddle 


-1 


unstable 


(iv) 


saddle 


-1 


unstable 


(v) 


saddle 


-1 


unstable 



rier is impossible in the traditional scalar field models 
pT| . The importance of the torsion cosmology is further 
promoted by this impossibility. In Fig. [2l we plot the dy- 
namical evolution of the equation of state Wef / for differ- 
ent initial values (iJ, <I>, R). Contrary to the quintessence 
and phantom model , the effective equation of state 



parameter Wef / is dependent on time that can cross the 
cosmological constant divide wa = — 1 from w^ff > — 1 
to Weff < — 1 as the observations mildly indicate. 

Critical points are always exact constant solutions in 
the context of autonomous dynamical systems. These 
points are often the extreme points of the orbits and 
therefore describe the asymptotic behavior. If the so- 
lutions interpolate between critical points they can be 
divided into a heteroclinic orbit and a homoclinic or- 
bit (a closed loop). The heteroclinic orbit connects two 
different critical points and homoclinic orbit is an orbit 
connecting a critical point to itself. In the dynamical 
analysis of cosmology, the heteroclinic orbit is more in- 
teresting [l^ . If the numerical calculation is associated 
with the critical points, then we will find all kinds of 
heteroclinic orbits. Especially, the heteroclinic orbit is 
shown in FiglSl which connects the positive and negative 
attractors. 

TABLE II: The physical properties of critical points for 
-1/9 < ai < 



critical points 


property 


™«// 


stability 


(i) 


saddle 


— CO 


unstable 


(ii) 


positive attractor 


-1 


stable 


(iii) 


negative attractor 


-1 


unstable 


(iv) 


saddle 


-1 


unstable 


(v) 


saddle 


-1 


unstable 
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FIG. 2: The evolution of the equation-of-state parameters 
w^ff for different initial values(Ji", $,ii) with —1/9 < ai < 0. 
We take ai = -1/10,6 = 8. 



V. CONCLUSION AND DISCUSSION 

In this paper, we investigate the dynamics of a torsion 
cosmology, in which we consider only the " scalar torsion" 
mode. This mode has certain distinctive and interesting 
qualities. We show that the late-time asymptotic behav- 
ior docs not always correspond to an oscillating aspect. 
In fact, only in the focus case can we declare that "scalar 
torsion" mode can contribute a quasinormal oscillating 
aspect to the expansion rate of the universe. There are 
only exact periodic solutions for the linearized system, 
which just correspond to the critical line (line of centers). 
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Via numerical calculation of the coupled nonlinear equa- 
tions Nester et al. @ plot that quasi-periodic solution 
near the focus. 

The late-time de Sitter attractor indicates that torsion 
cosmology is an elegant scheme and the scalar torsion 
mode is an interesting geometric quantity for physics. 
We show that the late-time de Sitter behaviors cover a 
wide range of the parameters and thus alleviate the fine- 
tuning problem. Furthermore, the torsion cosmology has 
considered the possibility that the dynamics scalar tor- 
sion (geometric field) could fully account for the accel- 
erated universe, which is naturally expected from space- 
time gauge theory. 



FIG. 3: The phase diagrams of {H, $, R) with -1/9 < ai < 0. 
The heterochnic orbit connects the critical points case (iii) to 
case (ii). We take ai = —1/10, 6 = 4. 
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